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ABSTRACT 

We present upper limits on the amplitude of r-mode oscillations and gravitational-radiation-induced spin-down rates 
in low-mass X-ray binary neutron stars, under the assumption that the quiescent neutron star luminosity is powered 
by dissipation from a steady-state r-mode. For masses <2 M Q we find dimensionless r-mode amplitudes in the 
range from about 1 x 1(Y 8 * * * * to 1 . * 5 x 1 0 6 . For the accreting millisecond X-ray pulsar sources with known quiescent 
spin-down rates, these limits suggest that <1% of the observed rate can be due to an unstable r-mode. Interestingly, 
the source with the highest amplitude limit, NGC 6440, could have an r-mode spin-down rate comparable to the 
observed, quiescent rate for SAX J1808— 3658. Thus, quiescent spin-down measurements for this source would be 
particularly interesting. For all sources considered here, our amplitude limits suggest that gravitational wave signals 
are likely too weak for detection with Advanced LIGO. Our highest mass model (2.21 M G ) can support enhanced, 
direct Urea neutrino emission in the core and thus can have higher r-mode amplitudes. Indeed, the inferred /--mode 
spin-down rates at these higher amplitudes are inconsistent with the observed spin-down rates for some of the 
sources, such as IGR J00291+5934 and XTE J 1751 — 305. In the absence of other significant sources of internal 
heat, these results could be used to place an upper limit on the masses of these sources if they were made of hadronic 
matter, or alternatively it could be used to probe the existence of exotic matter in them if their masses were known. 
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1. INTRODUCTION 

The /--modes are large-scale currents in neutron stars (NSs) 

that couple to gravitational radiation and remove energy and 

angular momentum from the star in the form of gravitational 

waves (GWs; Andersson 1998; Friedman & Morsink 1998; 
Friedman & Schutz 1978; Lindblom et al. 1998). The physics 

of these oscillations is important in relating the microscopic 
properties of dense matter — such as its viscosity and neutrino 

emissivity — to the macroscopic, observable properties of NSs 

such as their spin frequency, temperature, mass, and radius. 
The /--modes are unstable to gravitational radiation and their 

amplitudes can grow exponentially if viscous and other possible 
damping mechanisms are not large enough. Because they can 
affect the dynamic properties of the star — its spin frequency and 
temperature evolution — they are potentially important probes of 
the phases of dense matter. For example, the boundary of the 

r-mode instability region in the spin frequency-temperature 
plane, and in particular its minimum, which may determine 
the final rotation frequency of the star, is very different for stars 
with different interior compositions (Alford et al. 2012a; Ho 
et al. 2011; Haskell et al. 2012). 

Whether or not the r-mode instability limits the spin rates of 

some NSs, it is nevertheless important to explore its potential 
astrophysically observable signatures. While mass-radius mea- 

surements of NSs are important for constraining the equation 
of state (EOS) of dense matter, observations of their dynamic 
properties such as spin and thermal evolution are important and 
potentially more efficient in discriminating between different 
phases of dense matter. That is because dynamic properties are 
affected by the transport and thermodynamic properties of dense 
matter inside the star, such as viscosity, heat conductivity, and 


neutrino emissivity, which depend on low energy degrees of 
freedom and are very different depending on the phase of dense 
matter present (Alford et al. 2012a, 2012b; Alford & Schwenzer 
2012 ). 

Pulse timing observations of accreting millisecond X-ray 
pulsars (AMXPs) made with NASA’s Rossi X-ray Timing 
Explorer ( RXTE ) are beginning to reveal the long-term spin 
evolution of low-mass X-ray binary (LMXB) NSs. The pulsar 
recycling hypothesis, whereby millisecond radio pulsars acquire 
their fast spins via accretion, requires that at least some of these 
stars are spun up to hundreds of Hz — the current NS spin record 
being 716 Hz (Hessels et al. 2006). Still, it remains somewhat 
puzzling that the spin frequency distribution of AMXPs appears 
to cut off well below the mass-shedding limit of essentially all 
realistic NS EOSs (Chakrabarty et al. 2003), and it has been 
suggested that the r-mode instability may play a role in limiting 
NS spin rates (Andersson et al. 1999). Alternatively, recent work 
supports the idea that at least some of the AMXP population 
are in or close to “spin equilibrium,” essentially determined 
by the physics of magnetized accretion, and that an additional 
mechanism, such as the r-mode torque, is not required to halt 
their spin-up (White & Zhang 1997; Patruno et al. 2012b). 

While pulse timing noise in AMXPs has made the interpreta- 
tion of the X-ray pulsation data difficult, there are now several 
convincing measurements of spin-down during the quiescent 
phases between accretion outbursts. For example, both SAX 
J1808.4— 3658 (hereafter SAX J1808) and IGR J00291+5934 
(hereafter IGR J00291) show spin-down of the NS between out- 
bursts, at rates, v, of approximately 1-3 x 10 -15 Hz s -1 (Patruno 
2010; Hartman et al. 2009). The most convincing evidence for 
an accretion-induced spin-up during outbursts is in IGR J0029 1 , 
for which a peak value of about 3 x 10~ 13 Hz s _ 1 has now been 
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inferred (Patruno 2010). The magnitude of the “instantaneous” 
spin-up torque is larger in this source than the spin-down torque 
by about a factor of 100. The long-term evolution is then a 
competition between the accretion outbursts that spin it up and 
the quiescent spin-down intervals. If the outbursts are frequent 
enough then the star will be spun up; if not, then long-term 
spin-down occurs. As the spin-up torque is proportional to the 
mass accretion rate, the outcome can also be expressed in terms 
of the long-term, average mass accretion rate. 

The observed quiescent spin-down rate, v sl j, puts an upper 
limit on the torque that can be present due to any /--mode, as 
it has to be less than the observed torque of 2nlv s &, where 
/ is the moment of inertia of the NS. Now, spin-down of 
pulsars is typically ascribed to the magnetic-dipole torque that 
is almost certainly present at some level due to the large scale 
magnetic field of the NS. Indeed, magnetic field estimates are 
usually obtained by equating the entire observed spin-down 
to that expected theoretically for magnetic-dipole radiation, 
but the observed spin-down is due to the sum total of any 
torques present. Assuming that both /--mode and magnetic- 
dipole torques are present, then the magnitude of the r-mode 
torque is equal to the observed torque minus the magnetic- 
dipole torque. Since the magnetic field strengths of AMXPs are 
typically not known independently from spin-down estimates, 
the precise value of the magnetic-dipole torque is not known a 
priori. However, to the extent that the magnetic torque accounts 
for the majority of the observed spin-down, as is typically 
assumed, then the / -mode torque is likely to be much less than 
the measured spin-down torque. 

In addition to torquing the NS the viscous damping of the 
r-modes acts as an internal source of heat. Calculations of the 
coupled thermal and spin evolution of accreting NSs including 
the effects of r-mode heating and gravitational radiation have 
been carried out by several authors (Brown & Ushomirsky 
2000; Haskell & Patruno 2011; Ho et al. 2011; Haskell et al. 
2012). The primary assumption in these analyses has been that 
the long-term accretion (spin-up) torque is balanced by the 
r-mode torque due to the emission of gravitational radiation 
that carries away angular momentum, that is, the long-term 
average accretion torque is in equilibrium with the r-mode 
torque. As noted above, the recent pulse timing observations of 
AMXPs indicate that this equilibrium assumption is likely not 
realized in practice. Moreover, Brown & Ushomirsky (2000) 
showed that if the average fiducial accretion torque given by 
A acc = (M)(GMR) 1 ' 2 , where (M) is the long-term average 
mass accretion rate, were balanced by the r-mode torque, then 
the quiescent luminosities of some accreting NS transients 
should be substantially larger than observed due to heating 
from the r-modes. Here we use a similar argument to place 
upper bounds on the r-mode amplitudes that can be present in 
accreting NSs assuming that /--mode heating provides the source 
of NS luminosity in the absence of accretion. We then use these 
amplitude limits to assess the level of r-mode spin-down that 
can be present, its relation to observed spin-down rates when 
available, and the expected strength of gravitational radiation. 

The paper is organized as follows. In Section 2 we review 
the basic theory of the r-mode instability and how it couples 
to the thermal and spin evolution of NSs. In Section 3 we 
describe our methods for constraining the /-mode amplitudes 
using the observed properties of quiescent, LMXB NSs, and we 
assess the implications for r-mode spin-down and the emission 
of gravitational radiation. We provide a summary and discussion 
of the implications of our findings in Section 4. 


2. SPIN-DOWN DUE TO UNSTABLE r-MODES 

The “flow pattern” of the r-modes is prograde in the inertial 
frame and retrograde in the rotating frame, which means that 
it moves in the same direction as the star’s rotation as seen by 
an observer at infinity, but in the opposite direction as seen by 
an observer at rest on the star. Any mode that is retrograde 
in the co-rotating frame and prograde in the inertial frame 
grows as a result of its emitting GWs. This is the well-known 
Chandrasekhar-Friedman-Schutz mechanism, and it means that 
gravitational radiation drives the mode rather than dampening 
it. Viscosity, on the other hand, tends to damp the r-mode by 
transferring angular momentum from the mode to the rigid 
(unperturbed) star. The total angular momentum of a perturbed 
star can be written as 


j = m-j c (l) 

where I is the moment of inertia and £2 is the angular rotation 
frequency of the star and J c is the canonical angular momentum 
of the mode, given by 


3 _ 2 , 

J c = — J Mir Lie/ 2 
2 


(2) 


where a is the amplitude of the r-mode and / is a dimensionless 
constant defined by (Owen et al. 1998) 


J 


I 

MR 4 



(3) 


where p is the run of density within the NS. The moment of 
inertia of the star, 7, can also be written as I — I MR 1 , where I 
is defined by 


7 = 


8tt 

3MR 2 



(4) 


To derive the equations for the dynamical evolution of the 
star, we use the following argument by Ho & Lai (2000). The 
canonical angular momentum of the mode increases through 
gravitational radiation and decreases by transferring angular 
momentum to the star through viscosity 


dJ c _ 2 

dt Tg c 


2 

— Jc 
Ty 


(5) 


where the viscous damping time r y is given by (1/ty) = 
(1/tj) + (1/r b) + • ■ • • Here, r s and r# refer to shear and bulk 
viscosities, respectively, and the ellipsis denotes other possible 
dissipative mechanisms, such as boundary layer effects (Wu 
et al. 2001; Bildsten & Ushomirsky 2000). 

The second evolution equation is obtained by writing the 
conservation of the total angular momentum J of the perturbed 
star, which says that the total angular momentum of the star 
decreases due to gravitational radiation and increases due to 
accretion 


dJ 

dt 


2 J c 

Mice, 

tg 


(6) 


where /V acc is the accretion torque. For a fiducial torque /V acc can 
be written as M(GMR) 1 ^ 2 , which assumes that each accreted 
particle transfers to the star an angular momentum equal to the 
Keplerian value at the stellar radius R (Brown & Ushomirsky 
2000). In the previous equations the quantities 1/r, , where i is 
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Table 1 

Parameters of the Neutron Star Models 


Neutron Star 

Shell 

R 

(km) 

Qa- 

(Hz) 

/ 

j 

S 

V 

Cy 

L 

a 

s 

V 

e 

NS 1.4 M g 

Core 

11.5 

6020 

0.283 

1.81 x 10" 2 

7.68 x 10“ 5 

1.31 x 10" 3 

2.36 x 10“ 2 

1.91 x 10“ 2 

5 

3 

6 

1 

8 

NS 2.0 Mq 

Core 

11.0 

7670 

0.300 

2.05 x 10" 2 

2.25 x 10 -4 

1.16 x 10" 3 

2.64 x 10“ 2 

1.69 x 10“ 2 

5 

3 

6 

1 

8 

NS 2.21 Mq 

m.U. core 

10.0 

9310 

0.295 

2.02 x 10“ 2 

5.05 x 10“ 4 

9.34 x 10" 4 

2.62 x 10" 2 

1.29 x i tr 2 

5 

3 

6 

1 

8 


d.U. core 






1.16 x 10" 8 


2.31 x 10“ 5 


4 


6 


Notes. Radius, Kepler frequency, and radial integral parameters that appear in the moment of inertia, angular momentum of the mode, dissipative powers due to shear 
viscosity (from leptonic interactions) and bulk viscosity (due to Urea processes), specific heat, and neutrino luminosity for different neutron star models considered 
in this work (Alford & Schwenzer 2012; Alford et al. 2012a). For the 2.21 Mq model the bulk viscosity and neutrino luminosity parameters are different in the inner 
core where direct Urea processes are allowed; therefore, these values are given separately in the last row. 


either G, B, or S for gravitational radiation, bulk viscosity or 
shear viscosity timescales, respectively, are given by. 


Pi 

' 2 E r 


(7) 


where E c is the canonical energy of the /--mode, l\, is the power 
radiated by GWs, and Pb and Ps are the powers dissipated due 
to bulk and shear viscosity, respectively (in natural units where 
c = Ti = ks = 1; Owen et al. 1998; Alford et al. 2012a), 

E c = -a 2 Q. 2 JMR 2 
2 


\2m+2 


Pc = 


3271 (m - 1 Y m {m + 2) 

((2m + 1 )! !) 2 (m + \ ) 2m+2 “ m 

16m 


J,iGM 2 R 2m+2 a 2 Q 2m+4 


Pb = ~ 


VmAy^a^T 8 


(2m + 3)(m + 1 ) 5 k 2 


A 4 

2l EW 


Equations (5) and (6) and substituting J c from Equation (2). 
The third equation — which describes the temperature 
evolution — can be obtained by noting that the temperature of 
the star decreases due to thermal emission from the surface and 
neutrino emission from the interior, which in an average mass 
hadronic star is dominated by modified Urea processes (in a 
massive star direct Urea processes can also occur in the core 
and should be included in the neutrino emissivity as well), and 
it increases due to the viscous dissipation of the r-mode energy, 
Py. This gives the following equations for the evolution of spin 
frequency, r-mode amplitude, and temperature: 


(8) 


dQ. 

dt 

Qar A acc 

-2Q + -^ 

Xy I 


(12a) 

(9) 

da 

a 

a o 

— (1-crQ)- 

Xy 

a 

N acc 

(12b) 


dt 

re 

20 

/ 

(10) 

dE 

dt 

dt 

+ 

X 

1 

^4 

1 

II 

Pv 1 

+H, 

(12c) 


(m - l)(2m + l)5 m A£f n /? 3 a 2 Q 2 

Ps = Jn Q • (ID 

Here we consider m = 2 r-modes, and Aqcd and A E w 
are characteristic strong and electroweak scales introduced to 
make V and S dimensionless. In our calculations we have used 
Aqcd = 1 GeV and A E w = 100 GeV. The dimensionless 
parameters V , S, I and 7 , which involve radial integration 
over the star, and S and a are given in Table 1 for three 
different NS models that we study in this paper. All three models 
here are made of non-superfluid, hadronic npe matter with the 
Akmal-Pandharipande-Ravenhall (APR) EOS (Akmal et al. 
1998), which generates a reasonable NS mass-radius relation 
that is consistent with observational constraints (labeled AP4 in 
Hebeler et al. 2013 and Lattimer 2012), but they have different 
masses (1.4 Mq, 2.0 M 0 , and 2.21 Mq) and radii (Alford et al. 
2012a). The two models with masses of 1.4 Mq and 2.0 Mq 
only allow modified Urea neutrino emission in the core, but 
the one with a mass of 2.21 M Q allows direct Urea neutrino 
emission in a core of radius 5.9 km. Direct Urea processes 
are very sensitive to the proton fraction of dense matter. The 
required proton fraction is roughly 14% in the case of the APR 
EOS, reached at relatively high density n ~ 5«o, where no is 
the nuclear saturation density; this could be different for other 
EOSs (Alford et al. 2012a). 

The evolution equations for the amplitude of the r-mode, 
a, and spin frequency of the star, Q, can be written by using 


where Q = (37/27), and the viscous dissipated power is 
Py = Ps + Pb + ■ • • , and again the ellipsis denotes other 
possible dissipative processes. In Equation (12c), //represents 
other heating mechanisms that might be present in the star but 
are not related to the r-mode dissipation, such as deep crustal 
heating due to nuclear reactions in the NS crust (Haensel & 
Zdunik 2003; Gupta et al. 2007). Since our goal in this paper 
is to obtain upper limits on r-mode amplitudes we can safely 
set H = 0. We discuss this further in the next section. Here, 
L v , L y , and Cy, which are the total neutrino luminosity, photon 
luminosity, and specific heat of the star, are given by (in natural 
units) 


L v = 


4tt/? 3 Aq CD Z 


A 4 

21 EW 


Ly = 47r R a T* s 


Cy — 4tT Aq^ d R 2 CyT V , 


(13) 

(14) 

(15) 


where T is the core temperature of the star and 7/ n is the surface 
temperature. The dimensionless parameters L v and Cy (that 
involve radial integration over the star; Alford & Schwenzer 
2012) and 6 and v are given in Table 1 for different stellar 
models. 

Figure 1 shows the r-mode instability window computed for 
a 1.4 M q NS using the APR EOS (the modification to the 
instability window is rather modest in the case of the two other 


3 


The Astrophysical Journal, 773:140 (lOpp), 2013 August 20 


Mahmoodifar & Strohmayer 



Figure 1 . The r-mode instability region for a 1.4 M© NS constructed with the 
APR EOS in the spin frequency vs. core temperature plane. Also shown are 
some of the LMXBs that have been considered in this work. The horizontal line 
extending rightward from the temperature symbols (the black squares) shows 
the difference between two models for relating the surface temperature to the 
core temperature (i.e., the difference between a fully or partially accreted 
envelope). The difference between the core temperatures in these two cases 
gets larger as the surface temperature increases, but even for the highest surface 
temperature considered in this work the difference is not large enough to change 
our results. 

(A color version of this figure is available in the online journal.) 

stellar models considered in this paper (Alford et al. 2012a)). In 
the region above the curve, the gravitational radiation timescale 
is smaller than the viscous damping timescales; therefore, 
r- modes are unstable and their amplitudes grow exponentially 
in this region. As can be seen in this figure, most of the LMXB’s 
are in the unstable region for normal hadronic stars where the 
damping is due to shear viscosity from leptonic scattering (i.e., 
electron-electron and electron-proton scatterings, which are 
the dominant contributions to the shear viscosity of normal 
hadronic matter in NS cores) and bulk viscosity due to Urea 
processes; therefore, there must be some non-linear mechanism 
that saturates the amplitude of the unstable r- modes at a finite 
value. Suprathermal bulk viscosity (Alford et al. 2010) is one 
of these non-linear mechanisms, but in this case because only 
the core of the star is considered and the effects of the crust 
are ignored, it can only saturate the r-mode at large amplitudes 
(a ~ 1; Alford etal. 20 12c). Magnetohydrodynamic coupling to 
the stellar magnetic field is another mechanism that can damp 
the r-mode instability, but it can only saturate the r-mode at 
large amplitudes (a > 0.01) in the presence of a magnetic field 
significantly larger than ~10 8 G, which is characteristic of the 
LMXB sources considered here (Rezzolla et al. 2000, 2001). 
Mode coupling can saturate the r-mode at smaller amplitudes 
(i a ~ 10 4 ), but those values are still very large compared 
to the upper limits we find in this work (Arras et al. 2003; 
Bondarescu et al. 2007, 2009). At this point it is not entirely 
clear which mechanism is actually responsible for saturating the 
r-mode amplitude (none of the saturation mechanisms proposed 
thus far can saturate r-modes at the low amplitudes we find 
here); however, in this paper our primary interest is to obtain 


upper bounds on r-mode amplitudes from observations of NS 
transients and this does not require a precise understanding of 
the saturation mechanism. Understanding the detailed physics of 
the saturation mechanism is an important issue, but it is beyond 
the scope of this paper. 

Here we note that if we were to consider the existence of 
exotic forms of matter, such as strange quark matter in the star, 
then the shape of the instability window would change due to 
their different shear and bulk viscosities, and it is possible that 
the LMXBs considered here might fall outside of the resulting 
instability window (Alford et al. 2012a; Haskell et al. 2012; 
Schwenzer 2012). The shape of the instability window may also 
be different due to crust boundary layer effects (Ho et al. 201 1), 
but even in that case for realistic boundary layer models most 
of the LMXB sources are in the unstable region. Therefore, 
we assume that these sources have unstable /-modes that are 
emitting GWs, but the amplitude of the r-mode is not growing 
exponentially, and there is a mechanism that can saturate the 
growth of the r-mode. In addition, we assume that all of these 
stars are made of normal hadronic matter (constructed with the 
APR EOS) and all of the sources that are in the unstable region 
for hadronic stars in Figure 1 are emitting GWs at a constant 
amplitude. 

Since we do not know which mechanism is actually respon- 
sible for saturating the /--mode amplitude, we simply assume 
that there is a nonlinear mechanism that saturates the mode. 
When a hits the saturation amplitude, the right-hand side of 
Equation (12b) becomes zero, which implies that 


where we have neglected the last term in Equation (12b) since 
it is much smaller than the other terms. 

Therefore, when the amplitude is saturated, in all of 
the evolution equations (1/ty) should be replaced by 
(1/t g )(1/(1 - a~Q )) and \P V \ by (P G /(1 - a~Q)). Since vis- 
cosity alone cannot stop the growth of the r-mode and an 
extra mechanism is needed to do that, therefore at saturation 
1 / r G (1/(1 — a 2 Q)) is larger than (1/ty) and the reheating term 
on the right-hand side of the temperature evolution equation 
(Equation (12c)) will be larger than the reheating due to bulk 
and shear viscosity (Alford et al. 2012b). 

3. CONSTRAINING THE /--MODE AMPLITUDE 
3.1. Constraints from “Spin Equilibrium” 

Here we compare two different methods for constraining 
the r-mode amplitude, a, from observations of LMXB NS 
transients. The first one, which gives larger values for a, is based 
on the spin equilibrium assumption (Brown & Ushomirsky 
2000; Ho et al. 2011; Haskell & Patruno 2011; Watts et al. 
2008) where we assume that in an outburst-quiescence cycle 
all the spin-up torque due to accretion during the outburst is 
balanced by the r-mode spin-down torque due to gravitational 
radiation in the whole cycle. This is similar to the prescription 
considered by previous authors ( Brown & U shomirsky 2000; Ho 
et al. 201 1), but rather than using a “fiducial” torque estimated 
from the long-term average M, we can now use the observed 
spin-up rates and outburst properties to directly constrain the 
torque. Therefore we have 

2J C 

2jtIi>A= — - (17) 

to 
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Table 2 

Constraints from Spin Equilibrium 


Source 

A = 't r 

i^su 

(Hz s" 1 ) 

^sp.eq 

(1.4 M 0 ) 

^sp.eq 

(2.0 M 0 ) 

a sp.eq 

(2.21 M 0 ) 

IGR J00291 

13 

1363 

5 x lO^ 13 

1.46 x 10~ 7 

1.22 x 10~ 7 

1.41 x 10~ 7 

SAX J1808— 3658 

40 

2x365 

2.5 x 10“ 14 

3.20 x 10~ 7 

2.66 x 10~ 7 

3.08 x 10~ 7 

XTEJ1814— 338 

50 

19x365 

1.5 x 10" 14 

2.12 x 10~ 7 

1.76 x 10~ 7 

2.04 x 10~ 7 


Notes. Duty cycle A, spin-up rate v su (Patruno & Watts 2012), and upper limits on r-mode amplitudes from the “spin-equilibrium” 
condition Q'sp.eq^ for three neutron star transients with assumed masses of 1.4 M©, 2.0 M©, and 2.21 M©. 


where v is the spin-up rate during outburst and A = ( t a /t r ) 
is the ratio of the outburst duration, t 0 , to the recurrence time, 
t r . We can estimate the left-hand side of Equation (17) from 
X-ray observations of LMXBs and since the right-hand side is 
a function of a through J c (see Equation (2)), we can determine 
a. In Table 2 we give the results for a for three sources 
for which there are now estimates of the spin-up rate due to 
accretion. The values of a are given for the three different NS 
models considered in this work, and they are all in the range of 
^1-3 x 10 -7 . 

At these amplitudes, the inferred r-mode spin-down rate 
would be competitive with the magnetic dipole spin-down rate 
which almost certainly exists in these LMXBs and which is 
quite likely the dominant spin-down mechanism. Moreover, 
the amplitudes are also comparable to those deduced assuming 
r-mode spin equilibrium with the “fiducial” accretion torques 
estimated by Brown & Ushomirsky (2000). Those authors also 
demonstrated that at such amplitudes some of these objects 
should have significantly higher quiescent luminosities due to 
r-mode reheating than observed. These results suggest that the 
r-mode torque does not balance the accretion-driven spin-up 
torque and that r-mode amplitude estimates based on the “spin 
equilibrium” assumption will overestimate the true amplitude. 

3.2. Constraints from “Thermal Equilibrium” 

Here, we use the same thermal equilibrium argument out- 
lined by Brown & Ushomirsky, but rather than estimating the 
quiescent luminosity using the r-mode amplitude deduced from 
“spin equilibrium,” we use observations of the quiescent lumi- 
nosity of LMXBs to directly constrain the amplitude of the 
r-mode. This works because in a steady-state, gravitational 
radiation pumps energy into the r-mode at a rate given by 
Wj = (1/3)Q J c = — 2E c /tg ■ This expression has the familiar 
relationship for a power dissipated by an applied torque, and 
in this case it is simply the /--mode torque due to gravitational 
radiation. In a thermal steady-state, all of this energy must be 
dissipated in the star. Some fraction of this heat will be lost from 
the star due to neutrino emission and the rest will be radiated 
at the surface. It should be mentioned that the thermal steady- 
state is not an assumption but a rigorous result when the mode 
is saturated, and in particular it is independent of the cooling 
mechanism (Alford & Schwenzer 2012). We further assume that 
all of the energy emitted from the star during quiescence is due 
to the r-mode dissipation inside the star. This is equivalent to set- 
ting 77 = 0 in Equation (12c). The resulting r-mode amplitude 
limits are upper bounds in the sense that the observed luminos- 
ity reflects the contribution from r-mode heating as well as any 
additional sources of heat that are present, such as, for example, 
due to accretion and the nuclear processing of accreted material, 
so-called deep crustal heating (Haensel & Zdunik 2003; Gupta 


et al. 2007). If any such sources of heat are present, then the ac- 
tual r-mode amplitude will be less than the upper bounds given 
here. For the sources that we study in this work, since we know 
the values of surface temperatures and quiescent luminosities 
from observations, we can estimate the core temperature and 
therefore determine the neutrino luminosities to estimate the 
total amount of heat deposited in the core of these systems by 
gravitational radiation. 

To compute the core temperatures, we use Equation (A8) 
in Potekhin et al. (1997), which relates the effective surface 
temperature of the star 7/ n to the internal temperature 7), , which 
is the temperature at a fiducial boundary at pb = 10 10 g cm~ 3 * * * 
for a fully accreted envelope and is valid at not-too-high 
temperatures (7/ ^ 10 8 K), 

(^) 4 =(^)0^r <i8> 

where g = GM /( R 2 f I — r g / R) is the surface gravity and 
r g = 2 GM/c 2 . Here we assume that the NS’s core is isothermal 
and since the thermal conductivity of the crust is high (Brown 
& Cumming 2009), we have 7’ core = Tb to good approximation. 
If we instead use Equation (A9) in Potekhin et al. (1997), 
which gives Tb for a partially accreted envelope, and a column 
depth of ( P/g ) = 10 9 g cnT 2 (Haskell et al. 2012), we 
get core temperatures slightly higher than those derived from 
Equation (18). The right-hand side of the error bars on the 
temperatures in Figure 1 shows this difference. It is really only 
relevant for a single source, 4U 1608—522, but even in this 
case it is less than a 50% increase, and the difference is always 
small enough that it does not qualitatively change our results 
in the remainder of the paper. To compute 7/ we have used the 
effective surface temperatures, T e g , given by Heinke et al. (2007, 
2009) and Tomsick et al. (2004). ’ Note that those temperatures 
are computed for a 1.4 M Q NS with a 10 km radius and we 
have used the appropriate redshifts to compute 7’ core for our NS 
models. 

Having surface and core temperatures for these sources, we 
can use Equations (13) and (14) to evaluate their neutrino and 
thermal luminosities for different stellar models. The values 
of the core temperatures, as well as the photon and neutrino 
luminosities for different sources for the 1.4 M Q and 2.21 M e 
NS models, are computed and given in Table 3. Note that in 
the case of 1.4 M Q NSs the neutrino luminosity is only due to 
the modified Urea reactions, but for the 2.21 M 0 NS it is due 
to both direct and modified Urea reactions. By comparing the 

3 It should be emphasized that the temperatures given in Table 2 of Heinke 

et al. (2007, 2009) are effective temperatures at the surface of the star and not 

redshifted surface temperatures seen by an observer at infinity. We note that 

they have been incorrectly assumed to be redshifted temperatures in Degenaar 

et al. (2012) and Haskell et al. (2012). 
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Table 3 

Estimates of NS Temperatures and Luminosities 


Source 

V s 

(Hz) 

Distance 

(kpc) 

&T ff 
(eV) 

Tore (K) 

(1.4 M 0 ) 

Tore (K) 
(2.21 Mq) 

L y (erg s 3 ) 
(1.4 Mq) 

L v (erg s *) 
(1.4 Mq) 

L y (e rgs *) 
(2.21 M 0 ) 

L v (erg s 1 ) 
(2.21 Mq) 

4U 1608-522 

620 

4.1 ±0.4 

170 

1.25 x 10 s 

1.34 x 10 8 

1.21 x 10 34 

2.16 x 10 32 

3.03 x 10 34 

2.20 x 10 39 

IGR J0029 1+5934 

599 

5 ± 1 

71 

2.96 x 10 7 

3.18 x 10 7 

3.65 x 10 32 

2.06 x 10 27 

9.27 x 10 32 

3.87 x 10 35 

MXB 1659-29 

556 

11.5 ± 1.5 

55 

1.96 x 10 7 

2.07 x 10 7 

1.35 x 10 32 

7.64 x 10 25 

3.29 x 10 32 

2.96 x 10 34 

Aql X-l 

550 

4.55 ± 1.35 

94 

4.70 x 10 7 

5.06 x 10 7 

1.12 x 10 33 

8.40 x 10 28 

2.87 x 10 33 

6.35 x 10 36 

KS 1731-260 

524 

7.2 ± 1.0 

70 

2.89 x 10 7 

3.12 x 10 7 

3.44 x 10 32 

1.70 x 10 27 

8.88 x 10 32 

3.47 x 10 35 

XTE J1751-305 

435 

9 ± 3 

<71 

2.96 x 10 7 

3.18 x 10 7 

3.65 x 10 32 

2.06 x 10 27 

9.27 x 10 32 

3.87 x 10 35 

SAX J1808— 3658 

401 

3.5 ±0.1 

<30 

7.23 x 10 6 

7.69 x 10 6 

1.21 x 10 31 

2.63 x 10 22 

2.99 x 10 31 

7.78 x 10 31 

XTE J1814-338 

314 

6.7 ± 2.9 

<69 

2.82 x 10 7 

3.01 x 10 7 

3.24 x 10 32 

1.39 x 10 27 

8.12 x 10 32 

2.78 x 10 35 

NGC 6440 

205 

8.5 ±0.4 

87 

4.11 x 10 7 

4.46 x 10 7 

8.10 x 10 32 

2.88 x 10 28 

2.11 x 10 33 

2.97 x 10 36 

XTE J1807-294 

191 

8.35 ±3.65 

<51 

1.72 x 10 7 

1.83 x 10 7 

9.84 x 10 31 

2.71 x 10 25 

2.46 x 10 32 

1.44 x 10 34 

XTE J0929-314 

185 

7.8 ±4.2 

<50 

1.66 x 10 7 

1.79 x 10 7 

9.06 x 10 31 

2.06 x 10 25 

2.31 x 10 32 

1.23 x 10 34 


Notes. Spin frequency, distance to the source (Watts et al. 2008), effective temperature at the surface of the star (Heinke et al. 2007, 2009; Tomsick et al. 2004), core 
temperature, photon luminosity at the surface of the star, and neutrino luminosity for both 1 .4 Mq and 2.21 M© NSs. Note that kT t ff given in this table is for a 1 .4 M© 
NS with a radius of 10 km, but in computing the core temperatures and luminosities for different neutron star models the appropriate redshifts have been used. We 
note that for the sake of brevity the core temperatures for the 2 Mq models are not included in the table; however, their values are all ~5% less than those for the 
2.21 Mq model. 


neutrino and thermal luminosities in Table 3 one can see that 
at temperatures relevant for the LMXBs the standard neutrino 
cooling from modified Urea reactions is negligible compared 
to the photon emission from the surface of the star. However, 
if the star is massive enough to enable direct Urea reactions 
in the core (such as for our 2.21 M 0 NS model) then the 
neutrino emission will dominate the cooling process for surface 
temperatures higher than about 34 eV (Brown et al. 1998). 

The thermal equilibrium condition can be written as Wd = 
L v + L y , where reheating due to r-mode dissipation is given by 
Wd = {—2 E c /xgr) and is a function of r-mode amplitude, a. 
Therefore, a can be written in terms of luminosities as 


5 x 3 4 f L y +L v \ 1/2 

“ ~~ 2 8 /MR 3 £2 4 V 2 ttG ) 


(19) 


where L y = 4 tt R 2 o 7’jj ( is the thermal photon luminosity at the 
surface of the star. Here, R and 7’ c n are the stellar radius and 
surface temperature, respectively, and the neutrino luminosity 
is given by 


, _ 4jr ^ 3A QCD^DU 

v ~ A 4 

1 *-cx\r 


4tt A 3 A ( 


qcd-Cmu § 


A 4 

2V ew 


( 20 ) 


where T is the core temperature, 7 ?du is the radius of the 
core where direct Urea neutrino emission is allowed, and L 
is a dimensionless parameter given in Table 1. The thermal 
equilibrium condition for an NS with standard neutrino cooling 
(1.4 M g and 2.0 M 0 NSs in this study) can be approximated as 
Wd — L y , since the neutrino cooling in this case is negligible 
compared to the surface photon luminosity. 

As can be seen in Figure 1 , out of the 1 1 sources considered in 
this paper all but 2 of them are likely to have unstable r- modes, 
meaning that they are above the r-mode instability curve. The 
two most slowly rotating sources, XTE J1807 and XTE J0929, 
are outside the instability region for our NS models, which 
means they likely can no longer spin-down due to gravitational 
radiation from an r-mode. 4 Therefore, we only evaluate the 

4 It has been shown by Alford et al. (2012a) that the boundary of the 

instability region is insensitive to the quantitative details of the microscopic 

interactions that induce viscous damping in a given phase of dense matter. 


upper bounds on the r-mode amplitude for those nine sources 
within the instability window. Using Equation (19), we have 
evaluated a for all of those sources using the three different NS 
models considered in this work. The values of a are given in 
Table 4. As can be seen for the 1 .4 M 0 and 2.0 M 0 NSs, where 
there is no enhanced neutrino emission, the values of a range 
from 1 .07 x HT 8 to 1 .54 x ltU 6 , where NGC 6440— which has 
the lowest spin frequency — has the highest r-mode amplitude, 
as expected. In the case of 2.21 M 0 NSs, the upper bounds on 
a are larger, since due to the direct Urea neutrino emission W,/ 
can be larger, and a in this case ranges from 3.30 x 1 0 8 to 
8.03 x 10 -5 . 4U 1608—522 has the highest temperature among 
the sources considered, which implies a very large neutrino 
luminosity for the 2.21 M 0 NS model and as a result a large 
r-mode amplitude. The large value of the r-mode amplitude 
in this source may be ruled out by measurement of its spin- 
down rate, as is explained further in the next section. The high 
temperature of 4U 1608—522 could be explained if it has a 
lower mass, but it should also be noted that this system has been 
accreting for a long time and both its high quiescent luminosity 
and surface temperature may be due to the long-term accretion 
and our assumption in ascribing all the quiescent luminosity 
to the heat that comes from inside of the star may not be a 
good estimate for this system, but here we are only interested in 
obtaining upper limits on / -mode amplitude. 

3.3. r-mode Spin-down 

To see whether or not these results are plausible and what 
fraction of the quiescent spin-down of these sources can be due 
to gravitational radiation from r-mode oscillations, we use our 
results for a from the “thermal equilibrium” condition and insert 
them into the right-hand side of Equation ( 1 2a) to determine v 
in the absence of accretion (N acc = 0). The derived spin-down 
values for different NS models are given in Table 4, and are 
shown graphically in Figure 2 (vertical green lines). Comparing 
these results with the observed spin-down rates, which exist for 
IGR J00291, XTE 1751-305, and SAX J1808 (red diamonds 
in Figure 2), we find that in the case of 1 .4 M 0 ( x symbols) 
and 2.0 M 0 (square symbols) NSs, where there is no fast 
neutrino cooling present in the star, the / -mode spin-down can 
only provide about 1% of the observed spin-down rate, which 
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Table 4 

Upper Bounds on /--mode Amplitudes and NS Spin-down Rates 


Source 

a th.eq 

(1.4 M e ) 

^th.eq 

(2.0 Mq) 

^th.eq 
(2.21 Mq) 

v (Hz s 1 ) 
(1.4 Mq) 

v (Hz s 1 ) 
(2.0 Mq) 

v (Hz s 1 ) 
(2.21 Mq) 

v S d (Hz s 3 ) 
Observation 

4U 1608-522 

7.15 x 10~ 8 

6.60 x 10“ 8 

2.61 x 10~ 5 

-1.44 x 10- 15 

-1.78 x 10- 15 

-2.08 x 10” 10 


IGR J0029 1+5934 

1.41 x 10~ 8 

1.32 x 10“ 8 

3.99 x 10~ 7 

-4.42 x 10- 17 

-5.59 x 10- 17 

-3.82 x 10~ 14 

-3 x 10~ 15 

MXB 1659-29 

1.16 x 10^ 8 

1.07 x 10" 8 

1.49 x 10~ 7 

-1.78 x 10- 17 

-2.18 x 10- 17 

-3.16 x ur 15 


Aql X-l 

3.49 x 10^ 8 

3.27 x 10~ 8 

2.26 x 10“ 6 

-1.49 x 10~ 16 

-1.89 x 10~ 16 

-6.74 x 10~ 13 


KS 1731-260 

2.35 x 10“ 8 

2.20 x 10“ 8 

6.44 x 10~ 7 

-4.81 x 10" 17 

-6.09 x 10- 17 

-3.90 x 10" 14 


XTE J1751— 305 

5.09 x 10~ 8 

4.76 x 10~ 8 

1.44 x 10" 6 

-6.13 x 10" 17 

-7.74 x 10- 17 

-5.29 x 10” 14 

-5.5 x 10“ 15 

SAX J1808— 3658 

1.28 x 10~ 8 

1.19 x 10“ 8 

3.30 x 10~ 8 

-2.19 x 10" 18 

-2.74 x 10~ 18 

-1.57 x 10“ 17 

-5.5 x 10 16 

XTE J1814-338 

1.76 x 10~ 7 

1.67 x Hr 7 

4.49 x 10" 6 

-7.49 x 10- 17 

-9.73 x 10- 17 

-5.26 x 10 -14 


NGC 6440 

1.54 x 10^ 6 

1.45 x 10~ 6 

8.03 x 10~ 5 

-2.90 x 10~ 16 

-3.71 x 10~ 16 

-8.50 x 10~ 13 



Notes. Upper bounds on the r-mode amplitude from the "thermal equilibrium" condition that are consistent with quiescent luminosity data are given for different 
neutron star models. The gravitational radiation induced spin-down rates due to unstable r-modes as well as the observed spin-down rate for some of the sources are 
also given (Patruno 2010; Patruno & Watts 2012; Patruno et al. 2012a). 



i/ (Hz) 


Figure 2. Limits on the spin-down rates due to an r-mode torque for nine LMXB 
systems and a range of NS masses are shown in the v vs. v plane. The v limits 
for the 1.4, 2.0, and 2.21 Mq NS models obtained from the r-mode amplitude 
limits derived from observations of quiescent luminosities and temperatures 
(see the discussion in Section 3) are marked by the cross, square, and circle 
symbols, respectively. The vertical green lines connecting the symbols show the 
full range of v for each source (labeled). Also shown are two pairs of parallel 
lines representing magnetic-dipole (solid) and r-mode (dashed) braking laws. 
Lines are drawn for two values of the magnetic field, 10 s G (lower), and 10 9 
G (upper), as well as two values of a, 5 x 1 0 3 (lower) and 5 x 1 0 4 (upper). 
For systems with measured, quiescent spin-down rates, these values are marked 
with the red diamond symbols. For additional context, millisecond pulsars from 
the ATNF pulsar database (Manchester et al. 2005) are shown and denoted by 
the black plus symbols. 

(A color version of this figure is available in the online journal.) 

means that other spin-down mechanisms such as magnetic- 
dipole radiation are responsible for spinning down a 1.4M 0 
or 2.0 M q hadronic star with no fast cooling process. For the 
two AMXPs with relatively “slow” spin frequencies, NGC 6440 
and XTE J18 14 — 338, we find r-mode spin-down limits that are 
more competitive with observed values. Indeed, our limit for 
NGC 6440 is comparable to the measured, quiescent spin-down 
rate for SAX J1808, and thus spin-down measurements for this 
source would be particularly interesting. 

On the other hand, the r-mode amplitudes we obtain for 
the 2.21 Mq (circle symbols) hadronic model with direct Urea 


neutrino emission are only consistent with observations for SAX 
J1808 (and perhaps MXB 1659), as the inferred spin-down rates 
are either less than the observed rate for the source — in the case 
of SAX J1808 — or similar to the other observed rates, as for 
MXB 1659. For the remaining sources considered here, the 
2.21 Mq limits are likely not consistent with the observations, 
since such large amplitudes imply very large r-mode spin-down 
rates, and in the case of IGR J00291 and XTE J1 8 14 they are in 
fact larger than the observed values. If the neutrino luminosity 
from these sources was indeed as large as estimated with our 
2.21 Mq model, then in thermal equilibrium there must be a heat 
source that can supply it. Since the spin-down measurements 
for these sources indicate that r-mode heating (for this model) 
would be insufficient, several possibilities remain. First, there 
could be some additional source of heat other than r-mode 
dissipation that supplies the needed energy. However, we note 
that it would need to supply a substantial luminosity, as the direct 
Urea neutrino emission for this model outshines the photon 
luminosity by more than an order of magnitude, and we are not 
aware of any simple mechanisms that could provide the required 
luminosity. Second, the actual mass of these systems could be 
less than that of the model in question (2.21 Mq). Indeed, if it 
could be demonstrated that r-mode dissipation were the only 
mechanism that could produce such a large luminosity, then 
an upper limit on the mass would follow, and the limit would 
be the mass for which the neutrino plus photon luminosity 
matched the r-mode heating produced when the amplitude is 
large enough to produce a spin-down rate equal to the observed 
quiescent rate or our theoretical value for a high-mass NS 
model, whichever is smaller. This would be a conservative limit 
in the sense that it is likely that the r-mode torque does not 
account for all of the observed spin-down. Finally, our model 
assumptions, for example, the EOS and core composition, could 
be incorrect, with one possibility being the existence of exotic 
matter, such as kaon or pion condensates, or quark matter in 
the core, which have smaller neutrino emissivities than nucleon 
direct Urea processes, or if the pairing gaps for 3 I\ neutrons 
and 'So protons were larger than current theoretical values 
(this will be explained in more detail in the next paragraph). 
Interestingly, if the masses of these systems were known then 
one of the possibilities outlined above is precluded and then 
the observationally derived r-mode limits become sensitive to 
properties of the core, either the presence of exotic matter or 
perhaps additional heating physics. In this sense, further spin- 
down measurements and, where possible, mass constraints could 
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provide interesting new insights into the physics of dense NS 
matter. 

Our theoretical treatment of neutrino emission processes in 
this study, namely modified and direct Urea, spans the plausible 
range between “slow” and “fast” neutrino cooling processes. 
Here we have considered NS models made of non-superfluid 
hadronic matter with the APR EOS, but more realistically 
it is likely that neutrons and protons will be in a superfluid 
phase inside NSs. Therefore, a natural question would be 
whether or not our conclusions will still hold in the presence of 
superfluidity. Considering the presence of superfluidity in these 
sources, assuming that they are still inside the unstable region for 
r-modes, could have two possible effects on their cooling. The 
first one is neutrino emission due to Cooper pair breaking and 
formation (PBF) just below the superfluid critical temperature 
(Page et al. 2004) and the second one is the suppression of direct 
and modified Urea neutrino cooling at temperatures below the 
critical temperature. Here we explain why our qualitative results, 
and our argument about setting upper bounds on the masses 
of these sources, are not changed by considering the effect of 
superfluidity on the neutrino cooling of these sources. In our 
low-mass NS models ( M < 2 M Q ) where there is no fast cooling 
mechanism in the core, the thermal luminosity is much larger 
than the modified Urea neutrino luminosity, in fact by more 
than five orders of magnitude in all of the sources considered 
here except 4U 1608. Therefore, even if the temperature of these 
sources were just below the critical temperature of superfluidity, 
neutrino emission due to PBF, which is only about 10 times 
stronger than modified Urea neutrino emission, would still be 
much smaller than the thermal emission. Therefore, suppression 
of neutrino emission or neutrino emission due to PBF is not 
important for our low-mass NS models. What about high-mass 
NSs where fast cooling can happen in the core? Since direct Urea 
neutrino emission is much stronger than neutrino emission due 
to PBF, if direct Urea is not suppressed by superfluidity, it will be 
the dominant cooling mechanism. Now the question is whether 
or not neutron and proton pairing can suppress direct Urea 
neutrino emission in the core. Current theoretical results for 
the pairing gaps in 3 1\ neutron superfluid (Schwenk & Friman 
2004; Dong et al. 2013) and 'So proton superconducting phases 
(see for example Page et al. 2004; Yakovlev & Pethick 2004), 
which are relevant in the core of NSs, suggest that both of these 
pairing gaps are vanishingly small in the inner core of NSs where 
the direct Urea process can operate. Therefore, superfluidity is 
unlikely to suppress fast cooling processes in the core of NSs, 
and thus neglecting the effect of superfluidity does not change 
our qualitative results, assuming that /--modes are still unstable 
in the presence of superfluidity in these sources. We also note 
that the density profile (density versus radius) for an NS made 
of hadronic matter with the APR EOS is almost flat at the center 
of the star, which means that as the mass of the star increases 
above 2 M G (above which direct Urea processes can operate 
in the core of an NS made of hadronic matter with the APR 
EOS), there will be a sizable region in the core where direct 
Urea processes may happen, which can make it easier to obtain 
an upper limit on the mass of these sources using spin-down 
measurements. 

With typical values of a few x 10~ 8 , our derived amplitude 
upper limits suggest that for many LMXB NSs the /--modes 
are likely not excited to sufficient amplitudes to substantially 
affect their spin evolution. This begs the question of whether 
or not unstable, steady-state r-modes actually exist in these 
NSs. One possibility is that additional damping mechanisms, 


such as those perhaps associated with crust effects, such as the 
viscous friction at the crust-core boundary due to the coupling 
between core r-modes and crustal torsional modes (Levin & 
Ushomirsky 2001), superfluid mutual friction (Ho et al. 2011) 
or the existence of exotic matter in the core of NSs (Alford et al. 
2012a; Schwenzer 2012), are at work and modify the instability 
window so as to render these NSs stable to r-mode excitation. An 
interesting related question is whether the existence of r-modes 
at the amplitudes estimated here can be inferred directly from 
observations. Figure 2 shows both r-mode (dashed parallel lines) 
and magnetic-dipole (solid parallel lines) spin-down laws. The 
r-mode spin-down braking index, n = 7, is steeper compared 
to that for magnetic spin-down ( n = 3), thus at high-enough 
spin frequencies (well above a kHz) one might expect that the 
r-mode torque would eventually become competitive with or 
dominate the magnetic dipole torque. However, as of yet there 
are no known NSs spinning fast enough for this effect to become 
dominant, and depending on the EOS the mass-shedding limit 
might be reached before the r-mode torque becomes competitive 
with the magnetic torque. 

As discussed above, quiescent spin-down measurements have 
been typically attributed to the magnetic torque. For a number 
of sources considered here our r-mode amplitude limits support 
this presumption. Any spin-down contribution from an r-mode 
torque would be more easily identifiable if the magnetic field 
strengths of these NSs were constrained independently from 
the magnetic spin-down estimate. Moreover, identifying an 
r-mode spin-down would, in principle, be simpler for those NSs 
with the lowest magnetic torques, and thus field strengths. At 
present, the lowest inferred dipolar magnetic field strengths are 
~6 x 10 7 G. At this level the magnetic spin-down is of the order 
of u ~ 5 x 10 17 Hz s -1 , which is comparable to our derived 
r-mode spin-down limits for a number of sources considered 
here, assuming that their NSs are <2 M Q . In this regard, 
quiescent spin-down measurements for more of the sources 
considered here, in particular the AMXPs XTE J 1 8 14 and 
NGC 6440, would be extremely valuable. 

3.4. Gravitational Wave Amplitudes 

r-modes in NSs are one of the possible mechanisms for 
GW emission and they can be observationally interesting in 
newborn NSs and perhaps accreting NSs in LMXBs (Owen 
2010). Continuous GW emission from r-modes is dominated by 
l = m = 2 current quadrupole emission (Lindblom et al. 1998). 
The GW amplitude ho (strain tensor amplitude) is related to the 
r-mode amplitude a by the following equation (Owen 2010, 
2009) 

/ 87T G 1 -i t ~ 

h 0 = J T-aw 3 r MR 3 J, (21) 

V 5 c s r 

where r is the distance to the source, M and R are the mass 
and radius of the NS, J is the dimensionless parameter defined 
by Equation (3), and co r is the frequency of the r-mode, which 
is related to the angular spin frequency of the star £2 (for the 
m — l = 2 r-mode) by the following equation 

(Or » ^£2. (22) 

Using the upper limits on the r-mode amplitude of NSs in 
LMXBs derived above, we can obtain upper limits on the 
amplitude of the GWs emitted from these sources due to unstable 
r-modes. For the sources considered in this work, upper limits 
on the GW strain amplitude ho for the 1.4 M G and 2.0 M 0 NS 
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models are in the range of 1.8 x 10” 29 to 4.9 x 10” 28 , which 
is below the anticipated detectability threshold of Advanced 
LIGO (Abbott et al. 2010; Watts et al. 2008). In the case of 
the 2.21 M q NS model, since the r-modc amplitudes are larger 
than those for the low-mass NSs, we get larger values of ho, but 
even in this case for most of the sources ho is still below the 
detectability threshold of Advanced LIGO. The highest values 
of ho for a massive star are obtained for NGC 6440 and Aql X-l 
with an amplitude of the order of 8.5 x 10” 27 and 4U 1608—522 
with an amplitude of 1 .59 x 10” 25 . However, it should be noted 
that the large r-mode amplitudes in these sources, which would 
cause very large spin-down rates, may eventually be ruled out by 
future spin-down measurements. In this context it is important to 
restore the X-ray timing capability that was lost when RXTE was 
decommissioned in 2012 January. Missions planned or currently 
in development which could provide such a capability include 
India’s ASTROS AT (Paul 2012), ESA’s Large Observatory for 
X-ray Timing (Feroci et al. 2012), and NASA’s Neutron Star 
Interior Composition Explorer (Gendreau et al. 2012). 

The upper limits on GW amplitudes discussed here are related 
to the GW emission due to unstable / -modes, but our results do 
not exclude the possibility of having larger GW amplitudes 
in LMXBs from other GW emission mechanisms such as NS 
mountains (Ushomirsky et al. 2000; Haskell et al. 2006). It is 
worth mentioning that indirect upper limits on GW amplitude 
can be obtained for sources with observed spin-down rates, v, 
by assuming that all of the observed spin-down is due to GW 
emission (Owen 2010), 


h 


sd 

0 


l [45G1P 

8c 3 p 


(23) 


where P = (2jt /Q) is the observed pulse period and | P/P | = 
\i>/v\. Using this equation for the three sources with measured 
v and different NS models (i.e., different masses and radii) we 
obtain h s 0 d values that range from 4. 14 x 10” 28 to 6.53 x 10” 28 . 

4. CONCLUSIONS 

In this paper we have presented upper limits on r-mode 
amplitudes in LMXB NSs using their observed quiescent 
luminosities, temperatures, and spin-down rates. We calculated 
results for NS models constructed with the APR EOS (normal 
hadronic matter) with masses of 1.4, 2 and 2.21 M 0 , where our 
highest-mass model (2.21 M 0 ) can support enhanced, direct 
Urea neutrino emission in the core. We have used two different 
methods to calculate r-mode amplitudes. The first is based on 
the assumption that in an outburst-quiescence cycle all the spin- 
up torque due to accretion during the outburst is balanced by 
the r-mode spin-down torque due to gravitational radiation. 
This method gives amplitudes in the range of «l-3 x 10” 7 
for the sources with measured spin-up rates. Since in reality 
there are other sources of spin-down such as magnetic-dipole 
radiation that may be the dominant spin-down source, we use 
another method for computing the r-mode amplitude that does 
not ascribe all of the spin-down of the star to gravitational 
radiation and therefore gives tighter bounds on the amplitudes. 
This second method is based on the assumption that in a thermal 
steady-state some fraction of the heat that is generated in the 
star due to r-mode dissipation will be lost from the star by 
neutrino emission and the rest will be radiated at the surface. 
This assumes that all of the heat emitted from the surface of 
the star during quiescence is due to the r-mode dissipation 


inside the star, and thus provides an upper bound on the 
r-mode amplitude. We have computed core temperatures as 
well as neutrino and thermal (photon) luminosities for LMXB 
sources using measurements of the quiescent luminosities and 
surface temperatures and showed that at temperatures relevant 
for LMXB NSs, when there is no enhanced cooling mechanism, 
the cooling of the star is dominated by photon emission from 
the surface (for 1.4 and 2.0 Mq NS models), but in a massive 
star where direct Urea neutrino emission is allowed, the cooling 
is dominated by neutrino emission (for 7' c n > 34 eV). For the 
lower-mass NS models (1.4 and 2 M 0 ), we find dimensionless 
r-mode amplitudes in the range from about 1 x 10” 8 to 
1.5 x 10” 6 . We note that none of the saturation mechanisms 
proposed so far can saturate r-modes at these low amplitudes. 
Alternatively, the enhanced dissipation that would result from 
the existence of exotic matter in NS interiors could shift the 
instability window such that the LMXBs are perhaps stable to 
r-mode excitation (Alford et al. 2012a; Schwenzer 2012). 

For the AMXP sources with known quiescent spin-down 
rates these limits suggest that <1% of the observed rate can 
be due to an unstable r-mode. Interestingly, the AMXP with 
the highest amplitude limit, NGC 6440, could have an r-mode 
spin-down rate comparable to the observed, quiescent rate for 
SAX J1808. Thus, quiescent spin-down measurements for this 
source would be particularly interesting. Having enhanced, 
direct Urea neutrino emission in the core of our highest-mass 
model (2.21 M 0 ) means that the dissipated heat in the star can 
be larger and therefore it can have higher r-mode amplitudes. 
Indeed, the inferred r-mode spin-down rates at these higher 
amplitudes are inconsistent with the observed spin-down rates 
for some of the LMXB sources, such as IGR J00291 and XTE 
J 1751 — 305. If r-mode dissipation were the only mechanism 
available to produce this high luminosity, then this could be 
used to put an upper limit on the masses of these sources if they 
were made of hadronic matter. Alternatively, it could be used 
to probe the existence of exotic matter in them if the NS mass 
in these systems were known. In this way, future spin-down 
and NS mass measurements for the LMXB systems considered 
here, as well as for yet to be discovered systems, could open a 
new window on dense matter in NS interiors. For this as well 
as other reasons, we regard the re-establishment of a sensitive 
X-ray timing capability as vital to the use of NSs as natural 
laboratories for the study of dense matter. Using the results 
for / -mode amplitudes, the upper limits on GW amplitude due 
to /--modes have been computed. The upper limits on the GW 
strain amplitude ho for the 1.4 M 0 and 2.0 Mq NS models are 
in the range of 1.8 x 10” 29 to 4.9 x 10” 28 , which is below the 
anticipated detectability threshold of Advanced LIGO. In the 
case of the 2.21 Mq NS model, we obtain larger values for ho, 
but even in this case for most of the sources considered in this 
work, ho is still below the detectability threshold of Advanced 
LIGO. GWs due to other mechanisms such as NS mountains 
may have larger amplitudes in these systems. 
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